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Abstract. The spherical mean field approximation of a spin-1 model with p-body 
quenched disordered interaction is investigated. Depending on temperature and 
chemical potential the system is found in a paramagnetic or in a glassy phase and 
the transition between these phases can be of different nature. In given conditions 
inverse freezing occurs. As p — 2 the glassy phase is replica symmetric and the 
transition is always continuous in the phase diagram. For p > 2 the exact solution 
for the glassy phase is obtained by the one step replica symmetry breaking Ansatz. 
Different scenarios arise for both the dynamic and the thermodynamic transitions. 
These include (i) the usual random first order transition (Kauzmann-like) preceded by 
a dynamic transition, typical of mean-field glasses, (ii) a thermodynamic first order 
transition with phase coexistence and latent heat and (iii) a regime of inversion of static 
and dynamic transition lines. In the latter case a thermodynamic stable glassy phase, 
with zero configurational entropy, is dynamically accessible from the paramagnetic 
phase. Crossover between different transition regimes are analyzed by means of Replica 
Symmetry Breaking theory and a detailed study of the complexity and of the stability 
of the static solution is performed throughout the space of external thermodynamic 
parameters. 
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1. Introduction 

In macromolecular compounds in solution, complex molecules or polymeric chains can 
fold into practically inactive conformations, displaying a negligible interaction with the 
surrounding system. The presence of inactive components can induce the existence of 
a fluid phase at a temperature below the temperature range at which the system is 
in a solid phase (crystalline, semi-crystalline or amorphous, depending on the degree 
of frustration) [TJ El |3]. This corresponds to the occurrence of an inverse transition, 
else said "melting upon cooling", that is a reversible transition between a completely 
disordered isotropic fluid phase and a solid phase whose entropic content is greater than 
the entropy of the fluid. 

This effect can be reproduced and studied in statistical mechanical models on a 
lattice with bosonic spin-1 variables, where the holes s = play the role of inactive 
states. In these models the fluid phase is the paramagnet and the solid phase is either 
a ferromagnet (no or weak disorder) or a spin-glass (strong disorder). A prototype 
model with two-body interactions between bosonic spins is the Blume-Emery-Griffiths 
(BEG) model with ferromagnetic interaction [H [5} El [7] and its random magnetic 
implementation P El UTil HU H21 H31 HH UHl [16l HZl HEl [191 ED]. 

In presence of quenched disorder the random BEG model is known to display both 
a continuous paramagnet /spin- glass phase transition and a first order one. First order 
in the thermodynamic sense, i.e., with latent heat and a region of phase coexistence. 
Furthermore, inverse freezing takes place, with a spin-glass at high T and paramagnet 
at low T. These properties have been observed in the mean-field approximation [20J, 
where the self-consistent solution for the spin-glass phase is computed in the full Replica 
Symmetry Breaking (RSB) Parisi Ansatz [211 E21 E3], and on the cubic 3D lattice 
with nearest-neighbor couplings in systems exchanging two-body interaction, as well 
[2^1 [25] . The frustrated BEG model has been studied, as well, by means of numerical 
renormalization group techniques with apparently different results: on the Migdal- 
Kadanoff hierarchical lattice in 3D, it displays no inversion in the transition, nor a low 
T discontinuity |26j, whereas on the Wheatstone-bridge hierarchical lattice in dimension 
log 12/ log 2 ~ 3.585 the inversion seems to be there [27J. 

Mean-field spin-glass models with more than two-spin interactions, called p-spin 
models, are known to yield the so-called random first order transition, i.e., a phase 
transition across which no internal energy discontinuity occurs but the order parameter 
(the Edwards- Anderson overlap <7ea) jumps from zero to a finite value. Their glassy 
phase is described by an Ansatz with one RSB. The thermodynamic transition is 
preceded (in a cooling procedure) by a dynamic transition due to the onset of a very 
large number of metastable states separated by high barriers. The phenomenology of 
the |?-spin spin-glass systems is, in many respects, very similar to the one of structural 
glasses [281 [35] . These models are, therefore, sometimes called mean-field glasses. 

"Very large" means that the number of states M grows exponentially with the 
size N of the system: M ~ exp UN, where the coefficient E is the configurational 
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entropy, else called complexity in the framework of quenched disordered systems. This 
is a fundamental property both in mean-field systems [2HI ISH [3QI EH [32] and out of 
the range of validity of mean-field theory, e.g., in computer glass models [33], 151] . or, 
indirectly, by measuring the excess entropy of glasses in experiments, see, e.g., [35J and 
references there in. 

"High barriers" means that the free energy difference between a local minimum 
in the free energy functional of the configurational space (else called free energy 
landscape) and a nearby maximum (or saddle) grows with N. That is, it diverges 
in the thermodynamic limit. This is a specific artifact of mean-field glasses. The thus 
induced dynamic transition corresponds to the transition predicted by another mean- 
field theory for the dynamics of supercooled liquids: the mode coupling theory [36] . The 
thermodynamic transition occurring at a lower temperature is, instead, the mean-field 
equivalent of the so-called Kauzmann transition in glasses, else known as ideal glass 
transition, predicted by Gibbs and Di Marzio [37] . 

What happens to mean-field glasses when model features belonging to systems 
undergoing inverse freezing are introduced? The mean-field p-spin models are built 
either using discrete Ising variables [3_EJ [39] , soft spins [28] or spherical spins [30] , the 
latter two being approximation of Ising discrete variables that allow for an easier analytic 
treatment. 

In this work, we are interested in studying what happens when we mix the 
ingredients leading to an inverse transition (the hole state) and to a random first order 
transition (p-spin interaction), i.e., to provide a mean-field theory for inverse freezing 
between fluid and structural glass. We will see how, in this investigation, non-trivial 
features will arise, among which the inversion of dynamic and static transition and 
the consequent possibility of accessing low energy glassy states without running into 
dynamic arrest. 

2. Model 

The model Hamiltonian that we will consider, derives from 



with Si = 1,0,-1. The couplings J are Gaussian independent identically distributed 
variables with probability distributon: 



The coefficient D, known as crystal field in Blume Capel (BC) models [U [5j, plays the 
role of a chemical potential for the empty sites (holes). 

The p = 2 case was introduced in Ref. [8] and it has been studied subsequently 
throughout the years 0[IOl[IIl[I2[I3[Il[ig[II2 A preliminary study 

of the p-body case was presented in Ref. |40j. 




(1) 



ii<...<ip 




(2) 
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The spin-1 case can be translated into an Ising spin problem on lattice gas. Indeed, 
if we rewrite = o~iUi and crystal field D = D' + Tlog2, with o = ±1 and n = 0, 1 we 
have 

U = - ^ Ji 1 ...i p v il -..cTi p n il ...n ip + (D' -T log 2)^ ni (3) 

il<...<i p i 

We stress that the shift in chemical potential is necessary to keep the relative degeneracy 
of zero and non zero values of the spins in the partition sum [41 J. A different shift in D 
(or none) would, actually, define another model. Generally speaking a chemical potential 
transformation of the kind D = D' — T log 2r will provide the Ising spin on lattice gas 
model corresponding to a model with spin— 1 variables taking values s = 1 (or s = — 1) 
r times more frequently than value 0. E.g., r = 1/2 corresponds to s — 1, 0, 0, — 1; r = 1 
to the original case s — 1, 0, —1, cf. , Eq. (EJ); r = 4 to s — 1, 1, 0, —1, —1 and so on [6]. 
In the following we will address the variants with generic r, as well. 

In the Ising p-spin model, as p > 2, besides a random first order transition 
(RFOT) between a paramagnet and a spin-glass phase in the 1RSB Ansatz, also a 
lower temperature phase transition is expected to occur between the latter and a full 
RSB spin-glass phase. Indeed, this is what is known to occur in the D — > — oo limit 
[39J . Since we are exclusively interested in the transition between fluid and (mean- field) 
glass, thus 1RSB stable, we can simplify the computation by approximating the discrete 
spins with continuous real variables satisfying a global spherical constraint. 

As suggested in Ref. [42J, we introduce, to this aim, the variable r = a(2n— 1) = ±1, 
such that n = (err + l)/2, and eventually obtain the model Hamiltonian: 

1 N + 1 

U = ~Vv E Jn..,M,+^)---(o-> p + r tp ) + (D-T\og2)Y J ^ r - (4) 



2 p 

il<...<i v 



with 



N N 



5>* 2 = 5> 2 = iV - (5) 

i=l i=l 

Applying RSB theory j2U 1221 E3] we are going to investigate thermodynamics and 
complexity of the disordered model represented in Eq. (jl]). 

3. Replicated free energy and order parameters 

The replicated free energy, averaged over the distribution of disordered couplings, reads: 
U = --^ log JvQVTVR exp{-nNf3G(f3,D;Q,T,R)} (6) 



-nPG = - ((3D- log 2) £ + £ (Q ab + T ab + 2R ab f + 



ah 



1/1 \ n 

+ - In det ( - (QT + TQ) - R 2 \ + - log n + n, 
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l,n l,ri l,n 

a<b a<6 a<o 

where we have inserted the three overlap matrices 

NQ ab = e NT <* = E r " r *> ^ = E ^ ( ? ) 

i j i 

and, consequently, integrated over the spin variables. The saddle point, self-consistency, 
equations are: 

Qab = T a b = Rab, a b (8) 

-\[{QT- R>) - 1 R] ^ = ^^5 ab - ^ (Q ab + R^f- 1 0) 

From the free energy, cf. Eq. ([6]), in the zero replica limit 

/ = lim G(0, D, Q sp , T sp , R sp ) (10) 

one derives all thermodynamic quantities, such as the density of filled-in sites 

d = |- = limiy^±i; (11) 
dD n^on^ 2 ' V ; 

a 

the internal energy 

dpf r ^ J .. 2/3 1 



dp Dd -}^W^7X {Q ° t + T ^ + 2R ^ (12) 



and the entropy 



S = d log 2 + i log 7T + 1 



+ lim — 

«->o 2n 



2 " E (Gab + T ab + 2i? a6 ) p + lndet (QT - R 2 ) 



4(p+i) 



(13) 



In order to compute the above expressions the precise shape of matrices Q, T and 
i? has to be identified. There is no a priori method to deduce the correct form and one 
has to resort to an Ansatz. 

The simplest one is the Replica Symmetric (RS) Ansatz, where the discrete 
symmetry group S n , of permutation between replicas, holds for all two-indeces 
quantities. This means that the elements of each overlap matrix M only take one 
value y outside the diagonal and M can be parametrized as: 

M ab = (M d - M )5 ab + M , (14) 

where 5 is the Kronecker delta. Replica symmetry always holds for one index parameters, 
because physical properties of a single replica must be identical for each replica. Indeed, 
all single index quantities, like the diagonal part of the R overlap matrix, are index 
independent. 

RS is not always self-consistent, though. Studying the fluctuations in the space of 
replica matrices around a RS solution for Eqs. ((Hi)-© one finds that, for p > 2, the glassy 
solution is not stable and one has to break the symmetry between replicas in order to 
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find a self-consistent solution. To perform such stability analysis one needs the Hessian 



matrix. For our model it is computed in detail in Appendix A We will later report 
about the stability analysis of both the p = 2 case and the p > 2. Breaking the Replica 
Symmetry means to allow for different values in the off-diagonal elements of the overlap 
matrices. The symmetry can be broken step by step allowing a further overlap value at 
each step and organizing replicas in clusters with a hierarchical structure [2TJ [22] . For 
our purpose one step turns out to be sufficient, leading to a one step Replica Symmetry 
Breaking (1RSB) solution, in which a generic matrix takes the form: 

M ab = (M d - M x )5 ab + (M 1 - M )e ab + M , (15) 

where e ab = 1 if a and b belong to the same cluster of size m and e ab = otherwise. 

4. Thermodynamics for p = 2. Inverse freezing. 

For p = 2 the free energy evaluated in the RS Ansatz reads: 

- 2/3 lim G RS = -2(/3D - log 2)d + ^r{d 2 - q 2 ) + - + log(^) + log(4vr) + 2 (16) 

n— >0 2 Tj 

where 9 = 1 — d and rj = d — q are eigenvalues of the matrix QT — R 2 = Q 2 — R 2 . 
Self-consistency equations read: 

A = \ (17) 

/3 2 77=^ + 2(^-log2) (18) 
The RS solution turns out to be marginally stable with respect to fluctuations in the 



space of replica parameters, as shown in Appendix B In particular the lowest relevant 
eigenvalues of the Hessian, the so-called replicon, cf. Eq. (IB. 131) . is = 0. 

The self-consistency equation (|T7|) admits q = as a solution. This leads to a 
paramagnetic phase, characterized by a density dpM, determined imposing equation 
(flSflj ] and with a free-energy: 

- D) — — {pD - \og2)d PM + ^d 2 PM + 

+ l - log[(l - d PM )d PM ] + l - log(47r) + 1 (19) 

For low T and D the paramagnetic phase turns out to be unstable, in particular Eq. 
(IB. 131) becomes negative and a new solution to Eqs. (|T7|) - (|T8|) occurs with q > 0: 

d(T, D) = 1 - — — (20) 

v ' ; 2(l-D + Tlog2) v ; 

q(T, D) = 1 — T 21 

yv ' ; l-D + Tlog2 v ; 

| Actually, there is a region in the D — T diagram with three solutions for dpM, but only one turns 
out to be stable. See also Eq. P0)) with p = 2. 
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Figure 1. T, D phase diagram of the p = 2 model case. Inset: the same diagram 
is shown for models with different values or the parameter r = 1,2,4,8, ratio of 
the number of non-zero to zero values for the spin in the discrete counterpart of our 
spherical model. As r increases the inverse transition region is enhanced. 
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Figure 2. Density behavior at T = 0.4 across the SG/PM transition for the p = 2 
model. The paramagnetic density is shown (as a thin curve) also in its unstable phase. 



This SG solution is stable in the RS Ansatz in the whole region of the phase diagram 
where it exists, delimited by the transition line: 

D c (T) = l + Tlog2- 2(i ^ T) . (22) 

In Fig. [TJ we plot the phase diagram for the model. At the transition the overlap 
parameter q grows continuously from 0, the paramagnetic and the spin-glass solutions 
coincide and the density crosses continuously the transition, as shown in Fig. |2j This 
defines a second order transition, without latent heat or discontinuities in first order 
derivatives. As shown in the figure, a reentrance of the transition line points out the 
presence of inverse freezing. 

As we mentioned above, in Sec. El one can easily switch model using Schupper- 
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Figure 3. T, D phase diagram: in the inner region between at least one couple of 
complex eigenvalues is present in the Hessian in the space of replica fluctuations, as 
computed on the RS solution. In the other regions all eigenvalues are real. 



Shnerb variables [HI [7] controlling the degeneracy ratio r between filled and empty sites 
in the original model with discrete variables. The effect of increasing r is to increase 
the breadth of the interval in chemical potential for which inverse freezing takes place. 
In the inset of Fig. Q]we show the phase diagrams for some choices of r. 

Finally, we stress that, besides , cf. Eq. ( 1B.13j) . there are six more distinct 
eigenvalues of the stability Hessian, four of which are strictly real and positive. The other 
two, due to the n — > limit, can be complex conjugated and develop an imaginary part. 



Their real part is always positive. We refer to Appendix B for a detailed discussion. 
Here we only stress that, since the RS solution is exact for p = 2, the onset of imaginary 
eigenvalues is not a signature of lack of consistency of the replica Ansatz but an artifact 
due to the limit n — > in the replica calculation. In Figure [3] we plot the regions of the 
phase diagram in which imaginary stability eigenvalues occur. 

4-1. Random matrix approach 

To conclude the analysis of the p = 2 model we mention that, as shown in Refs. 
[1U H3j Hi], the thermodynamics can be computed also with the method applied by 
Kosterlitz, Thouless and Jones [15] to the spherical Sherrington-Kirkpatrick model, 
that is the D — > — oo limit of the present model. The method consists in describing the 
model in terms of the variables diagonalizing the interaction matrix: 

Jij^ = Jx^ (23) 

i 

The partition function turns out to be: 

Z = j dz 7T^ e ^{§(2*+f-riog2)-log|-^EA»««[AiWA 3 (*)]}_ ^ 
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where: 



A U z ) = o o ( 25 ) 



2 V 2 2 
» , x / £> -Tlog2\ 

A 2 (*) = £ U + ;r- iL - • (26) 



2 V 2 

In the thermodynamic limit, the sum over the eigenvalues can be evaluated as an 
integration having as measure Wigner's semicircle law: 

jr E — > f dJ P( J )f( J ) = f dJ^-V^Jf(J) (27) 

\ J — 2 J — 2 

Our replica calculation agrees with the result of Caiazzo et al. (121 H3j 03] . 
5. Static results for p > 2 

The situation is far richer, and interesting, when more than two-body interactions among 
spin- variables are considered. As p > 2, indeed, not only is the RS solution inconsistent 



in the SG phase (cf. stability analysis in Appendix B ) but the exact one step RSB 



solution yields different kinds of transitions, some of which atypical in the framework of 
mean-field glasses. For p > 2 the 1RSB free energy reads: 

_ pglRSB = _ (pp _ 1()g + ^ (d p + (m _ + (2g) 

If m — 1 1 

+ - 21og2 + log0 H log 770 H log 771 

2 \ m m 

where we have used the expression of the eigenvalues of the (Q 2 — R 2 ) matrix: 

9 = 1 - d (29) 

rj = d — q 1 (30) 

r}\ = d+ (m — - mg (31) 

We set g = 0, because of the absence of an external magnetic field. Self-consistency 
equations read: 

<K<h) = — (32) 

VoVi 

m - 0(gi) = + 2((3D - log 2) (33) 

z(y) = -. (34) 

p 



where 



and 



r>B 2 

<P(Q) = PJ Yq P ~ l (35) 

y ^ (36) 
m 

/ \ 1 - y + logj/ , . 

Z ^ = (l_ y )2 ( 3? ) 
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is the Crisanti-Sommers function [30] • The analysis of the stability of the 1RSB solution 



is reported in Sec. 15.61 and in |Appendix C 



The complexity functional, defined as the Legendre transform of the free energy 
functional, evaluated on the saddle point solution yielded by equations (152)1 and ( 133|) . 
logarithmically counts the number of metastable states [29]. It reads: 

E LT (m, D, P) = (3m 2 -^G(m, D, 0, qf(m, D, d s ?(m, D, (3)) (38) 

m 2 (3 2 (qf) p m^gf g r£ 
4 2^ p 2 g r)f 

To determine the dynamic transition line Td{D) one has to impose a maximal complexity 
(versus m) condition: 

^f = Vi = (p~l)Vo (39) 

in place Eq. (153]) . 

From the analysis of the paramagnetic solution (g = q± — 0), solving the self- 
consistency equation 

~d p (l-d) =T 2 (2d-l) + 2(T J D-T 2 log2)(l-d)d (40) 

one finds a region in the T, D phase diagram where three solutions for the PM density 
d occur. One of them is always unstable, whereas the other two coexist between the 
spinodal lines. The latter can be expressed, in a parametric form in d, by 

_ j~p~ d- 1+p / 2 (d + p-3dp + 2d 2 p + 2(1 - d) 2 d(p - 1) m(2r)) 
VP- 1 2(1 -d)y/4d?-4d + 1 

c/f/ 2 (l - d) 



T = VrtP^) /A ' (42) 

The first order transition line is, eventually, obtained by comparing the free energies of 
the two solutions dpM+ and dpM-- 

In Figs. H] and |5] we show the phase diagram for p = 3, both for r = 1 and 
r = 4. For r = 1 no reentrance of the thermodynamic transition line appears, as 
shown in Fig. HI but, increasing r, inverse freezing is clearly reproduced, cf. Fig. |5j 
The first order PM~/PM + line is the thick-dotted line and spinodal lines are plotted 
as thin-dotted lines. Qualitatively distinct transition regimes are identified. We will 
dedicate a subsection to each one of them, starting from the most conventional one. As 
a reference we follow the phase diagram in Fig. |4] moving from small (and negative) to 
large (positive) values of the chemical potential of the holes. 



5.1. Random First Order Transition (RFOT) 

For D small enough our model displays the same behaviour as the spherical p-spin 
model [30] and is an example of the well know RFOT. At high temperature the system 
is in a ergodic paramagnetic phase. As it is cooled down, crossing a temperature T d it 
undergoes a dynamical arrest, despite the static order parameter q\ is still null. Only at 
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Figure 4. T,D phase diagram of the spherical p — 3-spin BC model with filled-in to 
empty sites ratio r = 1. 
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Figure 5. T,D phase diagram of the spherical p = 3-spin BC model with filled-in 
to empty sites ratio r = 4. A reentrance occurs in the first order thermodynamic 
transition line. 



a lower temperature T s < a thermodynamic phase transition occurs. Critical static 
and dynamic lines are determined solving the equation systems ( ]32j) . (J33]), ( )34l) or ( l32l) . 
f )33|) . fl39|) . respectively. For p — 3, r = 1, e.g., in the D,T phase diagram they read 

„ /mx 1.8794T 1 / 3 - 2.49052T - 0.989769T 5 / 3 . . 

D '< r » = 1 - 02793^73 («) 

21.6337T 1 / 3 - 27.6822T - 11.5344T 5 / 3 
d( } ~ 12. - 16.6407T 2 / 3 ( } 

These curves are the m = 1-lines for statics and dynamics, respectively represented 
in Figs. HI |5] by the thick full (dark grey) and dashed (light grey) lines. As one can 
observe in Fig. El crossing the dynamical temperature from high T values, a collection 
of high complexity SG states arise at high free energy values. Each one of these states 
has a free-energy strictly greater than the free energy of the equilibrium state, namely 
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Figure 6. £(/; T, D = 0) in a region of temperature around static random first order 
(T S (D = 0) = 0.4739) and dynamic (T d (D = 0) = 0.4892) transition. 
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Figure 7. Detail of the D, T phase diagram with r = 1 with all significant critical 
points. A is the critical point of the PM + /PM~ transition; B is the merging point 
of the first order PM+/PM~ with the first order SG/PM^ line; E d>s is the limit of 
stability of the SG (dynamic/static) solution along the rn^.s = 1 line; / is the point of 
dynamic-static inversion of the m = 1 lines (where the f RSB solution is unstable) and 
^stab is the point of dynamic-static inversion in which the SG solution is stable; Fd, s 
are the points of dynamic-dynamic and static-static m-lines crossing, in the region of 
replica instability. 



the paramagnetic one. Indeed, from a strictly thermodynamic point of view, the stable 
state is the paramagnetic one, leaving each SG state as metastable. Decreasing further 
the temperature the system gains new states with lower free-energy and complexity, 
until, at T 
SG phase. 



T s , zero-complexity SG states appear, which turn out to yield the stable 
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5.2. Thermodynamic First Order Phase Transition 

At {Da, Ta) — (0.73, 0.44) a first order transition line starts developing in the PM phase, 
separating two distinct PM phases as T is lowered. The two phases are caracterized 
by their density, which has a jump crossing the transition line. We name PM± the 
paramagnetic phase at higher/lower density. At (Db,Tb) — (0.77,0.40) the first order 
PM + /PM~ transition line merges with a first order SG/PM" transition line D C (T), cf. 
Fig. [7J and the thermodynamic transition to the SG is not "random" anymore but, 
rather, a standard first order, ruled by the Clausius-Clapeyron equation, cf., e.g., Ref. 

The SG solution departs from the high-density PM + phase, i.e., SG and PM + 
belong to the same saddle point of Eq. fl29l) . see, e.g., Fig. [HJ For chemical potential 
values D < Db the thermodynamic transition to the SG is still the RFOT analyzed 
above. In Fig. [8] we plot the temperature behavior of the free energies of the PM~, 
PM + and SG phases in three qualitatively different cases: for D < Da, Da < D < Db 
and D > Db- 
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Figure 8. Free energy in the SG and PM phases in qualitatively different scenarios. 
The full curve represents $sg in the 1RSB Ansatz , ^pjyj-i- and ^p]yr- vs. T are 
plotted as dashed curves. D = 0.7,0.75,0.8,0.9. Top left: at D = 0.7 < D A , RFOT 
between PM - and SG, at T s . Top right: RFOT between PM+ and SG, at T s , and 
TFOT between PM" and PM+, at T c , for D = 0.75 (larger than D A and smaller than 
D B ). Bottom left: TFOT between PM - and SG, at T c , and RFOT between PM+ and 
metastable SG, at T s . Bottom right: RFOT between PM+ and metastable SG, at T s . 
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Figure 9. Left: Latent heat at the first order transition at D = 1 for the case r = 1. 
Inset: discontinuity of overlap q\ (dotted line) and density d (dashed line) across the 
same transition point at T c = 0.142. Right: Latent heat at D = 1.2 for the case 
r = 4, in presence of inverse freezing the SG phase transforms into the PM (at lower 
entropy As = sgc — s pm > 0) one as T decreases, taking heat from the thermal bath. 
Inset: the overlap q\ (dotted line) is different from zero above the critical temperature 
T c = 0.146. The density d (dashed line) decreases discontinuously in the low T PM 
phase. 



Along the D C (T) transition line, latent heat is taken from the glass in order to be 
transformed into a paramagnet, as we show in Fig. [9j 

Around that line the SG and the PM phases coexist and compete. Beyond the B 
point, the RFOT line T S (D) becomes the SG spinodal since on the right side of the 
D C (T) line the glassy phase is metastable with respect the PM one. That is, the whole 
hierarchy of global and local glassy minima lies above the global thermodynamically 
dominant PM minimum in the free energy landscape: &sg > &pm- 

This might seem weird in some respect. Indeed, in RSB theory, e.g., in the spherical 
p-spin model, the 1RSB SG solution stems out of the PM phase, which is RS. At a given 
T s a 1RSB solution appears with a discontinuous jump in q E \ and no discontinuity in 
the free energy As T < T s , the SG free energy is larger than the PM free energy but 
it is the one thermodynamically stable because of the n — 1, n - m or m - 1 factors 
present in the free energy expression Eq. ( )29|) whose sign is inverted in the n — y limit. 
This implies that looking for a minimal free energy at finite integer number of replicas 
corresponds to look for the maximal of the Parisi free energy. In our case, in fact, the 
SG 1RSB phase stems out of a PM (PM + ) phase, as shown in Fig. |HJ However, it 
competes with another phase PM~, corresponding to a different saddle point of Eqs. 
( l29l) . that does not transform into a SG and whose free energy does not involve overlap 
terms and, consequently, coefficients do not change sign in the replica calculation. 

As a homogeneous selecting rule, we can, thus, apply the minimal free-energy 
principle for all competing solutions with n > 1 and let n — > only after the choice of 
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Figure 10. Analysis of free energy in temperature at finite n = 0, 0.5, 1 , 1, 5, 2 in the 
RS approximation. The full curve is the SG solution in the zero replica limit. The 
dashed curves are finite n SG solutions. For n = 1 the SG free energy coincides with 
the PM + free energy. Each one of the SG curves crosses the free energy curve of the 
PM - phase signaling a first order phase transition between the PM~ phase, dominant 
at high T and the SG phase (dominant at low T). Since the SG free energy is computed 
in the RS Ansatz, the exact free energy is a little bit larger and the crossing point at 
slightly smaller T. 

the dominant stable phase has been done. To better exemplify this, we plot in Fig. [10] 
the RS approximated SG solution at finite n for different values of n and we compare 
it with both the free energies of PM~ and PM + . For n > 1 the free energy of the SG 
RS solution crosses the PM~ free energy for T lower than the transition temperature 
between the two PM phases, signaling the presence of a transition between the PM~ 
and the SG. The dominant phase is the one of least free energy, i.e., the SG phase. 
As n decreases the critical temperature decreases until, for n = 1, the SG free energy 
turns equal to the PM + one. As n < 1 the SG free energy shifts further towards small 
temperature but the thermodynamic order relationship with the free energy of the PM~ 
phase is unchanged. Breaking the replica symmetry on the same SG phase leads to a 
higher free energy and shifts the first order critical temperature towards smaller T. This 
does not affect the order relationship with respect to the other self-consistent solution, 
i.e., the PM~ phase. 

5.3. Dynamic- Static inversion 

As the chemical potential increases, the temperature interval between T S (D) and Td(D) 
decreases down to zero at T s (Dj) = T^Dj) = Tj = 0.3739, with Dj = 0.8905. For 
larger D the dynamic and static lines invert their position in T. Coming from high 
temperature, this means that the static line {m = 1 in Eq. ([3"1]) ) will be met before the 
dynamic one (m = 1 in Eq. (139]) ). According to the behavior of the complexity this 
implies that the lowest equilibrium glassy states are dynamically accessible at the static 
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Figure 11. £(/) for T = 0.35 at different values of D below the static transition. 



transition [§| and that excited glassy states only develop as T is decreased, as displayed 
in Fig. HU 

Actually, it can be seen from the complete stability analysis, cf . Appendix C that 
for D > DE d = 0.8737, slightly smaller than Dj, the SG dynamic solution spanned along 
the m = 1-line turns out to display a couple of complex conjugated eigenvalues with 
negative real part. This implies that for D > De the first thermodynamically stable SG 
solution occurs with the static parameter m < 1 and zero complexity for the excited 
states. The static-dynamic inversion then occurs, self-consistently, for values of m < 1, 
at Dj = 0.877, Tr h = 0.377. 



5.4- m-lines inversion, compressibility and replica stability 

As one can see from Fig. [121 solving equations ( 1321) . ( 1331) and Eq. ( 13^1) [else (1391)]. 
and following solutions at fixed m in the D,T plane, each T m (D) curve ("m-line") 
develops a maximum and T m (D) lines at different m values cross. This leads to an 
ambiguity: there is a region in the phase diagram in which in each point two different 
solution of Eqs. fl32|) . fl33]) .f l34|) occur. A criterium is, then, required to select the correct 
(metastable) phase in each point and, furthermore, a way to mark the spinodal of the 
SG. The solution comes from the positiveness of the compressibility k and from the 
complete analysis of the stability of the 1RSB solution in the replica space. Moreover, 
as we already mention, the m-line inversion comply with the fact that the spinodal SG 
lines (both dynamic and static) will cease to be a line of constant m = 1. 

§ The RFOT static transition line is actually a spinodal, since the thermodynamic dominant phase is 
PM beyond the D C (T) line. 
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Figure 12. T(D) lines of constant m (from static equations). As D is large enough 
they bend and cross each other. Left: complete behavior of m-lines. Right: only curves 
along which the compressibility is non-negative are plotted. An equivalent situation 
occurs for the dynamic m-lines. The thick continuous curve is the line where complex 
stability eigenvalues acquire negative real part: solutions on the right hand side of the 
stability line are rejected. 



5. 5. Compressibility 

As a thermodynamic quantity the compressibility is usually defined as: 
<91ogV^ 



k 



dp 



(45) 



In our model, where the role of the pressure is played by the chemical potential D and 
the inverse specific volume is the density of filled-in states d, the compressibility reads: 



k = d 



dd~ 



dD 



1 dd 
"d dD 



T 



1 dd_ 

d dm 



T 



dm 



(46) 



where the last equality comes from the Fubini's implicit function theorem. Obviously, 
solutions with negative compressibility are unphysical and have to be discarded. At 
fixed T the value m = m yielding the zero point of the compressibility is given by 



/ dm 




\dD 









(47) 



using, e.g., Eq. ( )33|) for the D(m) dependence. Beyond m(T) the compressibility turns 
negative and that solution will be rejected. This kind of "selection rule" has to be used 
for both the static and the dynamic m-lines. 

The positiveness of k is certainly a necessary condition for thermodynamic stability. 
On top of that we anticipate that, from the replica stability analysis, we found a small 
region of the phase diagram where k > but the real part of a couple of complex- 
conjugated eigenvalues among A 2) i, or among A31 of the Hessian in the replica space for 
the 1RSB Ansatz, cf. Eqs. flC.19j) . (IC.24j) . is negative (see also the right panel of Fig. 
T2i) . The spinodal lines plotted in Figs. HJ |5] are drawn taking the latter observation 
into account. See also the right panel of Fig. 
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Figure 13. Compressibility vs, T at D = 0.5 < Da, in the RFOT regime (left) and 
at D = 1 > D B in the TFOT regime (right). 



The previous discussion implies that the dynamic SG spinodal is the m = 1-line 
on the left side of the point F d , (D Fd = 0.98, T Fd = 0.33) and a line of decreasing m on 
the right side. Also the static SG spinodal undergo the same change at the point F s 
(Dp s = 1.03, Tp g = 0.32). We remark, however, that this change in behavior arises for 
D > Dj stah , where the m = 1-line cease to be stable in the replica space and D > D t , 
after that the static-dynamic inversion occurs. 



5.6. Stability of 1RSB solution 

In order to determine the above phase diagram it has been necessary to study the 



stability of the 1RSB solution in the T, D plane. In Appendix C we show the details 



of our analysis. The phenomenology looks similar to the RS case with both real and 
complex eigenvalues. The stability of the SG solution is ruled by {, the lowest real 
eigenvalue of the matrix Hn, cf. Eq. (IC.llj) . Evaluated on the static SG solution it is 
always positive, thus, according to our analysis, the static SG 1RSB solution is always 
stable. If, however, the static Eq. ( 13"1|) is not imposed and the parameter m is, thus, 
left undetermined, the positiveness of A^ 3 ] is equivalent to the condition: 

r]i > {p - 1)770 (48) 

and the equality A^ 3 ] = turns out to coincide with the dynamical condition Eq. ( )39|) . 
The other real eigenvalues are positive in the whole phase diagram. 

In the region of phase coexistence, near the SG spinodal, the stability analysis 
actually revelas more subtle features. Indeed, as in the p = 2 case, cf. Sec. HI 
complex eigenvalues are present and some of them develop a negative real part. In 
particular, solutions with negative compressibility always yield an imaginary eigenvalues 
with negative real part. Moreover, in a small region of the T,D plane solutions with 
n > also display complex eigenvalues, A 2 ,i or A 3j i, cf. Eqs. ( 1C.19I) . ( 1C.24I) with 
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Figure 14. Region of complex eigenvalues in the T, D plane for the p — 3 model. 



negative real part. As anticipated, in determining the phase diagram we have thus 
rejected solutions with negative replica eigenvalues, including all those with negative 
compressibility. 

6. Conclusions 

We have been analyzing, at the level of mean-field theory, a family of models mixing the 
properties of (i) Blume-Capel models with quenched disorder and (ii) spin-glass models 
with many-body (p > 2) interactions. The first ones are known to reproduce the inverse 
freezing occurring, e.g., in complex macromolecular compounds in solution [HE} [3], and 
display first order thermodynamic phase transitions even in presence of strong disorder, 
The latter yield most of the basic features of the standard folklore of glassy systems and 
are often refereed to as mean-field glasses. Our aim has been to investigate to behavior 
of dynamic arrest, displayed in mean-field glasses, as a system undergoes an inverse 
thermodynamic transition. 

We have been studying the problem by means of replica theory for the two 
qualitatively distinct cases with p = 2 and p > 2, both looking at the thermodynamic 
and at the dynamic properties. The dynamics has been tackled by means of the analysis 
of the complexity of the free energy landscape. The model approximation used is to 
pass from discrete spin-1 variables to spherical ones (i.e., continuous variables with an 
overall "spherical" constraint) keeping the original probabilistic relationship between 
states of the discrete spin: r = [# of times s = 1]/ [# of times s = 0]. This allows for 
thoroughly analytical computation of complicated physical quantities. We have been 
mainly focusing on the case r = 1, i.e., the spherical formulation of the original BC 
model but we have considered, as well, model cases with different r values. Indeed, in 
those cases where the spherical analogue of the original spin-1 system (r = 1) do not 
show inverse freezing, a simple variant of the model for larger r is suitable to describe 
such phenomenon. 
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The external parameters driving possible transitions are the temperature and a 
chemical potential for the empty states of the spin-1 variables, the so-called "crystal 
field" D. For D = — oo one has no holes and the limit of the model is the spherical 
spin-glass. In that limit, in the p = 2 case [15] the spin-glass phase is described by a 
replica symmetric Ansatz and the transition turns out to be always continuous, whereas, 
as p > 2 the RS solution is unstable and the right physics is described by a one step 
Replica Symmetry Breaking phase [30]. On top of it the transition is not continuous 
(at zero external magnetic field) in the order parameter and it is often refereed to as 
random first order. 

In the present study we find that, for p — 2, the situation is not very different for the 
original spherical spin-glass, with a spin-glass phase stable under the RS Ansatz. The 
only addition is the presence of a transition between competing paramagnetic phases at 
high T (D) and the occurrence of inverse freezing for large enough chemical potential 
D. For p > 2 the situation is, instead, drastically modified. Besides the standard 
RFOT, we find that, for large enough D, a first order thermodynamic phase transition 
between the spin-glass and the paramagnetic phase takes over, with latent heat and 
jumps in density and overlap order parameter. The spin-glass becomes metastable with 
respect to the paramagnet and the RFOT line plays the role of the spin-glass spinodal. 
Moreover, analyzing in detail the structure of the solutions in the phase diagram and 
their stability properties, we find that the spin-glass phase can be approached, in given 
regions, without incurring in dynamic arrest, unlike in the standard p-spin and that the 
first spin-glass solution display values of replica symmetry breaking parameter m less 
than one. The study of the complexity functional clarify how this observation amounts 
to a scenario in which, in a cooling experiment, lowest glassy minima of the free energy 
landscape develop first and excited metastable glassy states only arise as temperature 
is further lowered. This provides a model case where actual dynamics might be followed 
along given paths of the phase diagram into the deep glassy phase without being stuck 
in the threshold states. 
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Appendix A. Hessian of the spherical p-spin Blume-Capel model 

In order to verify the stability of a saddle-point calculation, the positiveness of the 
fluctuation matrix evaluated on the stationary solution has to be checked. This Hessian 
matrix is constructed as the second order variation of the free-energy potential ((6]) with 
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respect to the three overlap matrices ([7j) and, indeed, depends on two couples of indexes: 

»<G>2£=*£|jr. O,0 = Q,S,T, (A.1) 
In detail the elements read 



dG 2 1 

- B ab 5 a , c 5 bd + - [A aiC (TAT)bd + 2(TA) a>c (AT) bd + 



where: 



dQabdQ 



cd 



+(TAT) atC A bd 
dG 2 „ , „ I 



dQ ab dT cd 



- B ab 5 a J bd + - [A a ^TAQ) bd + (QA) ac (AT) bd + 
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+{AT) ac (AQ) bd + (TAQ) ac A bd \ - ~ (A 6 A C + A ac 5 fe 
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dT ab dQ cd 
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To check its positiveness of the Hessian, one needs to compute, and determine the 
sign of all the eigenvalues of its action on fluctuation matrices (<5o's): 



A 5o ab . = ( 

cd \ 



d 2 G 



dO ab dQ 



led + 



d 2 G 



cd 



B ab 5o ab 



0O ab dT cd 
IP 



<ttcd + 



d 2 G 



dO ab dR cd 
log LAI I 5o cd . 



8r 



cd 



(A.4) 



(A.5) 



In order to solve Eq. (1A.5I) . a set of 3n ~ n coupled equations, one needs to find a 
way to project the system on proper subspaces, that allows to translate the calculation 
to an usual eigenvalues problem. 



Appendix B. Stability of the RS solution 

In order to understand the appearance of complex eigenvalues, we start the stability 
analysis retaining n finite. In the RS ansatz the matix (1A.3j) becomes: 



A -1 



A9r] 



x 



ab 



x5 ab + w 



w = 



B^ b = -4>'(q) 



/3 2 p(p-l) . 



-Q 



p-2 



B 



B n 



B n 



(B.l) 
(B.2) 
(B.3) 
(B.4) 
(B.5) 



where 8 = 1 — d, rj = d — q and r\\ = d — (1 — n)q, with r/i — > rj if n — > 0. The parameter 
x is always positive. For later convenience, we define: 

D ++ = x 2 (t] + 9f- H ++ = xw[r] 2 + 46 2 + r](A6 - r]^}; (B.6) 

= (r] + 6)(r] - 6); = xw(t] 2 - 46 2 - 7]^); (B.7) 

= (rj - Of- H^ = xw{rf + A6 2 -r l {2e + r ll ) (B.8) 



Projection on subspace S\ 

Fluctuations that involve the overlap of one replica with the others, are selected by 



projecting the Hessian (lA.lj) in the subspace Si defined by: 

5r a , a = 0, 22 Sqa '-r = 



with dimension d(Si) = d to t — 4n 
the Hessian matrix reduces to: 

/ i 



3n(n— 3) 
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displaying no explicit n dependence. The characteristic polynomial, evaluated on the 
SG solution, reads: 

P[H X \ = (A - |) (A 2 - hA + d) 
3 _ {p _ 2)02) 



6i 



Cl 



p-2 

8r] 2 6 2 



Interaction with p = 2 

For p = 2ci = 0, 6i>0: two eigenvalues will be positive and one will be zero. 



(i) 



Ai 

aS 2) = 
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We have marginal stability of the RS-SG solution. 



(B.10) 

(B.H) 
(B.12) 



(B.13) 
(B.14) 

(B.15) 



Interaction with p > 2 and replica symmetry breaking 

If p > 2, instead, C\ < 0, leading to a negative zero of the polynomial and causes the 
instability of the RS-SG solution. As p > 2 one has to resort to a more complicated 
Ansatz, the 1RSB scheme of computation, whose stability will be analyzed in the 
following appendix. We now continue the analysis of the RS fluctuations in other 
subspaces. 



Projection on subspace S 2 

Fluctuations of the cluster of replicas are selected by projecting in S 2 , defined as: 

^5r 7i7 = 0, ^5g 7j(5 = 

7 7,<5 



with d(S 2 ) = 4n — 4. In S 2 the Hessian reduces to: 
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++ 



n— 2 tt 




where: 



with 



H[:- 



n—2 tt 

n—2 tt 
— H +- 

~\ H +- 



n—A TT 
~ H +- 4 



H__ + w 



A 



(4) 











A 4 


4B d 


+ x - 



H__ +w J 

(B.16) 



(B.17) 



(B.18) 



Random, thermodynamic and inverse first order transitions 24 

For generic n, apart from n = 4, H 2 is not symmetric, so the matrix can have complex 
eigenvalues. The characteristic polynomial of the H 2 matrix, evaluated on the SG 
solution, reads: 



P[H 2 ](A) = (-A+ | + ^j^^j [-A 3 + a 2 (n)A 2 - b 2 (n)A + c 2 ( 



n 



(B.19) 



In the limit n — > 0, the coefficients are: 



a 2 = —2wf] 2 x 3 

b 2 = -wx + 6x 4 r) 4 (B.20) 
c 2 = 5x 2 r] 2 — 2wx9 2 

Interaction with p = 2 

For the case p = 2, a numerical study ensures that all the solutions of the characteristic 
polynomial are real positive or complex conjugated with a positive real part. The region 
of the (D,T) plane where the polynomial has complex zeros (shown in Fig IB} can be 
determined imposing: 

- + — >0 (B.21) 
4 27 K ' 

where (in terms of the polynomial coefficients): 

" "\~i C2 (B.22) 

Stability of the paramagnetic phase 

In the paramagnetic phase w oc q = 0, implying H±± = 0, cf. Eqs. f ]B.6j) - f)B.8j) and the 
Hessian simplifies as H 2 = H[. Its characteristic polynomial is: 

P\H 2 ] = (A - A 4 ) (A - |) (-A 2 - b 2 A + c 2 ) (B.23) 

where now, cf. Eq. ( IB. 101) : 



3 / 2d' 1 ■ 2d ■ 1 )2 2 



^-lelli^W-^J (R24) 

c 2 = ^ 3^ (B-25) 



32(1 -d) 2 d 2 

Its zeros are: 



AS 2) = \ (B.26) 
A?> = ~ h + ^ = 4C2 (B.27) 



Af = ' U2 '^- 4C2 (B.28) 
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The first one is strictly positive for any p. 

The second one is always positive for p > 2 and ensures the stability of the PM 
solution in the whole phase-space. For p = 2, instead, Af M marks the PM solution as 
unstable in the region where a SG solution exists and permits to discard an unphysical 
solution preventing a PM-PM first order transition. A^ M > A^ M and, indeed, it is 
irrilevant. 

The fourth eigenvalue, equal to the derivative of Eq. (1331 with q — 0, allows to 
discard another unphysical solution among the PM solutions. 

Projection on subspace S3 

In the subspace S 3 = (Si U S 2 )' L (d(S 3 ) = 4) the Hessian reduces to a matrix 
if 3 = if 2 + 0(n); indeed there are no new eigenvalues. 



Appendix C. Stability of the IRS solution 

The paramagnetic solution is not affected by the Ansatz choice; its stablity analysis 



is indeed the same as the one in Appendix B Moreover for p = 2 is not necessary to 



break the replica symmetry, so in this section we choose p > 3 without loss of generality. 
These observations allow us to work directly with Eq.'s (I32II33I) imposed. Instead, in 
order to study also dynamical solution Eq. ( 1341) is not imposed. 
In the 1RSB ansatz the matix (1A.3I) becomes: 

A -l $ab Vl-VO _ r . (n , X 

x = J- ( C - 2 ) 
4% 

2m9rj rji ^ ^ 

(C.4) 

B a ^b = ^'(qi)e ab = ^(P-I) gr 2 £a6 s B ^ (C 5) 

?2„ 



B„ t = ^ Bd (C.6) 

o4 

where e a j, is the matrix introduced in section [3j For later convenience, is usuful to define: 

2 

D ±>± = ^(r ] o±e)(r ]0 ±e) (C.7) 
£± ± = — (7/o±^)(r7o ±0) + 



2 V^-"/V/u--, ^ 

_ x 2 (r]! — ?7o) 2 , wx(r)i - r]o){ M r 1 ^ + Vo + |0?i ~ Vo)) 

U±± = A 2 1 

4m z m 
+ — [( m ± e)( Vo ±6) + ((±1 ± 1)6 + 2 Vo )( Vl - 770) + (77! - Vo ) 2 ] 
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and some linear combinations: 



Et^ = D++ + mE. 



x(x + mw) 



Gq = D±± + 2mE±± + m 2 G±± 
± := 2E±± + mG±± 



( Vl + r ]o ±26)(r h +r ]o ±2e) 
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Projection on subspace Si 
Si,o 

In the subspace S^o, defined by 

Sr aa = 0, (e 5q) ab = (Sq t) ab = 0, Sq ab e ab = 
with dim(Sifi) = ^^^(n — m)(m 2 — 2m), the Hessian reduces to: 

D ++ D ++ -l -2D+. \ 



Hi, 



D. 



* Jl 

2 "I 



-2D, 



-2D, 



-2D + _ 
AD__ + x 



Eigenvalues of H 
i,o 2 



1.0 



Ai 



(2,3) 



2L>_ 



+ T ± 



4 



2L>_ 



D ++ + 



3x 



8Di 



are always definite positive. 



(C.8) 

(C.9) 
(CIO) 



Si ; x In the subspace Si t i, defined by: 

o~r aa = 0, (e 5q) ab = (Sq e) ab = 0, Sq ab (l - e ab ) = 
with dim(Si t i) = 3|(m — 3), 



( D A 



H ltl = 

the eigenvalues are: 

A?> = £ 

1,1 2 

Afi 3) = 2L>_ 
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B qi —2D + - — 2B qi 



\ 



-2D, 



2B qi -2D + _ - 2B qi 4D__ +x- AB gi J 



D, 



x 



35, 



91 



(C.ll) 

(C.12) 
(C.13) 



± 



2D^ ~D ++ + ^ - D, 



8p+_ - B, 



qi. 



(1 2 3) (1) 

Since Ai X o ; and \ are positive definite, the stability condition is expressed by 



A\{ > 0, which reduces to: 

7/i > (p - 1)7/0 

which concides with the dynamical condition, Eq. (I31?|) . 



(C.14) 
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Projection on subspace S 2 

$2,0 

In the subspace S-^o, orthogonal to Si,- U S^i defined by: 

5r aa = 0, (e 5q e) ab = 0, (5q e) ab e ab = 
with dim(S2fl) = 3^ (— — l) (m — 1), the Hessian reduces to: 



H- 



2,0 



jtH — h 2x+mw 
^0 4 



771++ 2x+mw 171+- 

4 — "^0 



V 



-2E n 



-2EI 

4S n ~ - 



/ 



The characteristic polynomial is: 



with: 



P[H 2 , ] = (a - (A» - 6 20 A + c 20 ) 

3 [4£ 2 + (1 + yfrft 



J2Q 



C20 



64y6 2 r)j 

(i + y) 2 
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(C.15) 



(C.16) 



128y 2 9 2 r]f 

with y = T]o/r]i, cf. Eq fl36|) . The positiveness of both 6 20 an d c 2 o implies positive zeros 
for the polynomial. 



In the subspace S^,-, orthogonal to Si )0 U SV,i and restricted by: 

(e5re) a6 = J^(5r cc e ac ) = 0, (e5qe) ab = 0, (5qe) ab (l - e ab ) 



with dim(S2,i) = 3—(m — 1) +n — ^. Defining the matrix: 
#2,1 := 



#1,1 








AD_. 


_ + x- AB qi 



(C.17) 



(C.18) 



the Hessian projection in S 2 ,i reduces to: 

/ E ++ E ++ 



H 2 ,i = H' 2 l + (m - 2) 



V 



F — 1 
-2E+- 



-2E- 



-2E 
-2E 

4E— 



-2E + - 
-2E+- 



m-2 



-4 



m-2 



f§ AE— + 1 

8E — +y 8E — +y 

m—2 m—2 



J 



(C.19) 



This matrix has to be studied numerically on the static (and dynamic) 1RSB solution. 
In the whole D, T plane it displays two real postive eigenvalues. The other two can 
be both real (and positive) or complex conjugated. In the region of D, T values where 
complex eigenvalues occur they have a positive real part, apart from a small subregion 
next to the static (resp. dynamic) RFOT transition line. 
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Projection on subspace S3 



3.0 



In the subspace S^o, orthogonal to all other subspaces and defined by: 

5r aa = 0, (e 5q e) ab e ab = 
with dim(Szfl) = 3^ — l), the Hessian reduces to: 



(C.20) 



H 



3,0 



C++ 

— h x+m 

<Jq 2 

-2G+- 



C++ x+mw 

C++ 

-2G+- 



-2G C 
-2G; 



\ 



4Gr, + 2x + raw 



The characteristic polynomial is: 
P[H 3fi ] = (a 

with: 

3/1 
030 = — I — + — 

c 30 



2x + mw\ . . 2 , . . 
J (A — 6 30 A + c 30 ) 



(C.21) 



(C.22) 



16 \Q 2 ' vl 
1 



329 2 r] 2 

the positiveness of both 6 30 and C30 implies positive zeros for the polynomial. 



S 3 ,i 

In the subspace S^i, orthogonal to all other subspaces and resricted by 

(e 5q e) ab (l - t ab ) = 
with dim(S 3t i) = 4— the Hessian reduces to: 

/ Gf + Gt + - 

C++ y c++ 

(jt x — g Lr 1 

-2Gf- -2G+- 



H 3A = H> + (m - 1) 



V 



4 G ^ _ 4^1 



-2G+- -2G+- 
-2Gt~ -2G+- 
4G~+y 4G~+3/ 



" m— 1 



m—1 



m—1 



m—1 



(C.23) 



(C.24) 



Also this matrix has to be studied numerically, as Eq. (1C.19j) and we evaluated its 
values point by point in the phase diagram. This matrix, as well, has two real positive 
definite eigenvalues together with two complex conjugated eigenvalues (in a region of 
the D,T plane). Some points, next to the RFOT transition have complex eigenvalues 
with negative real part. Comparing with the study of compressibility, cf. Sec. I5.5[ all 
point with a negative compressibility also have at least two complex eigenvalues with a 
negative real part. The opposite, as we have shown in. Fig [HI is not always true. 
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